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[^ ■ Abstract: We study in detail the power spectra of scalar and tensor perturbations 

generated during inflation in loop quantum cosmology (LQC). After clarifying in a 

t^^ ■ novel quantitative way how inverse-volume corrections arise in inhomogeneous set- 

■ tings, we show that they can generate large running spectral indices, which generally 

lead to an enhancement of power at large scales. We provide explicit formula for 
the scalar/tensor power spectra under the slow-roll approximation, by taking into 

^ ■ account corrections of order higher than the runnings. Via a standard analysis, 

j^ ■ we place observational bounds on the inverse-volume quantum correction 5 oc a~°" 

(a > 0, a is the scale factor) and the slow-roll parameter Cy for power-law potentials 
as well as exponential potentials by using the data of WMAP 7yr combined with 
other observations. We derive the constraints on 5 for two pivot wavenumbers k^ for 
several values of 5. The quadratic potential can be compatible with the data even 
in the presence of the LQC corrections, but the quartic potential is in tension with 
observations. We also flnd that the upper bounds on 5(/co) for given a and ko are 
insensitive to the choice of the inflaton potentials. 
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1. Introduction 

Both the construction of quantum gravity and the question of its observational tests 
are beset by a host of problems. On the one hand, quantum gravity, in whatever 
approach, must face many mathematical obstacles before it can be completed to a 
consistent theory. On the other hand, assuming that a consistent theory of quantum 
gravity does exist, dimensional arguments suggest that its observational implications 
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are of small importance. In the realm of cosmology, for instance, they are estimated 
to be of the tiny size of the Planck length divided by the Hubble distance. 

Between the two extremes of conceptual inconsistency and observational irrele- 
vance lies a window of opportunity in which quantum gravity is likely to fall. It is 
true that we do not yet know how to make quantum gravity fully consistent, and it is 
true that its effects for early-universe cosmology should be expected to be small. But 
in trying to make some quantum-gravity modified cosmo logical equations consistent, 
it has been found that there can be stronger effects than dimensional arguments sug- 
gest. Consistency requirements, especially of loop quantum gravity, lead to modified 
spacetime structures that depart from the usual continuum, implying unexpected 
effects. 

To partially bridge the gap between fundamental developments and loop quan- 
tum gravity phenomenology and observations, one considers effective dynamics com- 
ing from constraint functions evaluated on a particular background and on a large 
class of semiclassical states. In generally- covariant systems, the dynamics is fully 
constrained, and the constraint functionals on phase space generate gauge transfor- 
mations obeying an algebra that reveals the structure of spacetime deformations. The 
algebra of these gauge generators thus shows what underlying notion of spacetime 
covariance is realized, or whether covariance might be broken by quantum effects, 
making the theory inconsistent. If a consistent version with an unbroken (but per- 
haps deformed) gauge algebra exists, it can be evaluated for potential observational 
implications. 

These effective constraints and their algebra, at the present stage of develop- 
ments, can be evaluated for perturbative inhomogeneities around a cosmological 
background. While loop quantum gravity is background independent in the sense 
that no spacetime metric is assumed before the theory is quantized, a background and 
the associated perturbation theory can be introduced at the effective level. In the case 
of interest here, the background is cosmological, a flat Friedmann-Robert son- Walker 
(FRW) spacetime with a rolling scalar field and perturbed by linear inhomogeneous 
fluctuations of the metric-matter degrees of freedom. The idea is to implement as 
many quantum corrections as possible, and study how the inflationary dynamics is 
modifled. In this context, it is important that no gauge flxing be used before quan- 
tization, as such a step would invariably eliminate important consistency conditions 
by flat, not by solving them. The result would be a framework whose "predictions" 
depend on how the gauge was flxed in the flrst place. 

This program, challenging as it is, has been carried out only partially so far, and 
in gradual steps. First, the full set of constraints was derived for vector [|l|, tensor 
[^ , and scalar modes 0, ^ in the presence of small inverse- volume corrections. The 
gravitational wave spectra have been studied in P, 0, where the effect of inverse- 
volume corrections and their observability were discussed. The scalar inflationary 
spectra and the full set of linear-order cosmological observables was then derived in 



-2- 



[§] , thus making it possible to place observational bounds on the quantum corrections 
themselves p (see also refs. |10| for early related works). The reason why these 



studies concentrate on inverse- volume corrections is mainly technical; in fact, the 
closure of the constraint algebra has been verified only in this case (and only in the 
limit of small corrections). A class of consistent constraints with a closed algebra 
is known also for vector and tensor modes in the presence of holonomy corrections 
[0, ^ ^, and since recently also for the scalar sector, where anomaly cancellation is 
more difficult to work out |[Tl| . Inspections of cosmological holonomy effects have so 



far been limited to the tensor sector |g, |T^, |I3|, |14|, [L5 |. 

Exactly isotropic minisuperspace models, where the situation is reversed and 
holonomy corrections are easier to implement than inverse-volume corrections, pro- 
vide another reason why it is interesting to focus attention on inverse-volume cor- 
rections. Effective equations available for certain matter contents with a dominat- 



ing kinetic energy [|I6|, ^ suggest that holonomy corrections are significant only 



in regimes of near-Planckian densities |TB| but not during the timespan relevant 
for early-universe cosmology, including infiation. Inverse-volume corrections, on the 
other hand, do not directly react to the density but rather to the discreteness scale 
of quantum gravity, which is not determined immediately by the usual cosmological 
parameters. The question of whether they are small or can play a significant role 
must be answered by a self-consistent treatment. 

Such a treatment shows that inverse-volume corrections present an example of 
quantum-gravity effects that can be larger than what dimensional arguments sug- 
gest . Here we present the full details of the analysis briefly reported in for a 
quadratic inflationary potential, enriching it with new constraints on other poten- 
tials. From a cosmological perspective, we shall provide the complete set of slow-roll 
equations as functions of the potential, extend the likelihood analysis to quartic and 
exponential potentials, and discuss how the experimental pivot scale and cosmic 
variance affect the results. 

Before examining the details and experimental bounds of the model, from a 
quantum-gravity perspective we will clarify some conceptual issues which must be 
taken into account for a consistent treatment of inverse-volume corrections. In par- 
ticular, we justify for the flrst time why inverse- volume corrections depend only on 
triad variables, and not also on connections. Until now, this was regarded as a tech- 
nical assumption devoid of physical motivations. Here we show it as a consequence 
of general but precise semiclassical arguments. Further, we spell the reason why 
inverse-volume corrections are not suppressed at the inflationary density scale. 

The plan of the paper is as follows. We begin in section ^ by discussing inverse- 
volume corrections in LQC and their justiflcation in inhomogeneous models, provid- 
ing the flrst implementation for a class of semiclassical states sufficiently large to 
be used in effective equations. The lattice refinement picture has been introduced 
and developed in a number of papers, but a major twofold open issue remains. On 
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one hand, there is the need to justify why inverse-volume corrections depend only 
on triad variables and not on holononiies. On the other hand, it is not clear how 
the fundamental discreteness scale, previously introduced in media res, arises in the 
lattice refinement framework. Section ^ does both systematically for the first time. 
The relation between (the size of) inverse- volume and holonomy quantum corrections 
is clarified in section |2.3| , while section ^.4| is a recapitulation of past criticism and 
a discussion of how it is addressed by the present arguments. After that, we turn 



to an application for observational cosmology. In section ^^1 we review the formulae 
of the cosmological observables in the Hubble slow-roll tower [Q . In section |3]^ we 
reexpress these quantities in the slow-roll parameters as functions of the inflationary 
potential. In section ^ the observables are recast as functions of the momentum pivot 
scale for a ready use in numerical programs. The effect of cosmic variance on the 
scalar power spectrum is also discussed therein. In section ^ we shall carry out the 
likelihood analysis to constrain the inverse-volume corrections in the presence of sev- 
eral different infiaton potentials by using the observational data of Cosmic Microwave 
Background (CMB) combined with other datasets. 

2. Cosmology with a discrete scale 

One of the main features of loop quantum gravity (LQG), shared with other ap- 
proaches to quantum gravity, is the appearance of discrete spatial structures replac- 
ing the classical continuum of general relativity. It is often expected that the scale of 
the discreteness is determined by the Planck length ipi = \^Gh, but if discreteness 
is fundamental, its scale must be set by a dynamical parameter of some underlying 
state, just as the lattice spacing of a crystal is determined by the interaction of atoms. 
In this section, we develop the cosmological picture of dynamics of discrete space, 
highlighting the form of quantum corrections to be expected. Readers more inter- 
ested in potentially observable consequences may skip this technical part, but those 
acquainted with LQG will find a fresh discussion and justification of inverse-volume 
corrections and the lattice refinement picture. 

In loop quantum gravity, such states are represented by spin networks, graphs 
in an embedding space whose edges e are labeled by spin quantum numbers jg. 
The quantum number determines the area of an elementary plaquette intersecting 
only one edge e, given by ^ = '~fipi\^jeije + !)• As the plaquette is enlarged, its 
geometrical size changes only when it begins to intersect another edge, increasing in 
quantum jumps. In the area formula, 7 is the Barbero-Immirzi parameter, whose 
size (slightly less than one) can be inferred from computations of black-hole entropy 
[0, 1^. As expected, the scale is set by the Planck length for dimensional reasons, 
but the actual size is given by the spin quantum number. Its values in a specific 
physical situation have to be derived from the LQG dynamical equations, a task 
which remains extremely difficult. However, given the form in which j^ appears in 
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the dynamical equations, its implications for physics can be traced and parametrized 
in sufficiently general form so as to analyze effects phenomenologically. 

2.1 Scales 

In order to model this situation in general terms, we begin with a nearly isotropic 
spacetime and a chunk of space of some comoving size Vo, as measured by the ex- 
tensions in some set of coordinates. The geometrical size is then V = Voo^, where 
a is the scale factor. We complement this classical picture with a discrete quantum 
picture, in which the same chunk of space is made up from nearly-isotropic discrete 
building blocks, all of the same size v. If there are M discrete blocks in a region of 
size V, we have the relationship v = Vqo? /M . The elementary volume f , or the linear 
scale L = v^^^, will be our main parameter, related to the quantum state via its labels 
je- The elementary quantum-gravity scale L need not be exactly the Planck length, 
depending on what je are realized. Instead of using the je, which are subject to com- 
plicated dynamics, it turns out to be more useful to refer to L in phenomenological 
parametrizations. Similarly, we define the quantum-gravity density scale 

_ 3 

which equals S/Stt times the Planck density for L = ip\. 

In loop quantum gravity, the discreteness is mathematically seen as a rather 
direct consequence of the fact that the fundamental operators are holonomies along 
curves e, computed for a certain form of gravitational connection, the Ashtekar- 
Barbero connection A^,^ while the connection itself is not a well-defined operator. 
For a nearly isotropic spacetime, there is only one nontrivial connection component, 
given by c = 7(2 in terms of the proper-time derivative of the scale factor. The 
classical holonomies are h^ = V exp{J^d\e"'Al^Ti), with Tj = iaj/2 proportional to 
Pauli matrices and path ordering indicated by V. Every he takes values in the 
compact group SU(2), whose representations appear as the spin labels of edges je, 
giving rise to discrete conjugate variables. 

Another consequence of one being able to represent only holonomies, not con- 
nection components, is that the usual polynomial terms in connection-dependent 
Hamiltonians are replaced by the whole series obtained by expanding the exponential 
expression for an holonomy. In this way, higher-order corrections are implemented 
in the dynamics. Corrections become significant when the argument of holonomies, 
given by line integrals of Al^ along the spin-network edges, is of order one. For a 
nearly isotropic connection A^ = C(5^, the integral along straight lines reduces to EqC, 
where io is the coordinate (i.e., comoving) length of the edge. If the edge is elemen- 
tary and of the discreteness size of our underlying state, we have £o = L/a = v^^^ /a, 



^Indices a,b, ■ ■ ■ — 1, 2, 3 run over space directions, while i,j,--- — 1, 2, 3 are internal indices in 
the su(2) algebra. 
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and the condition for holonomy corrections becoming large is v^/^c/a ~ 1. More intu- 
itively, holonomy corrections become large when the Hubble scale H~^ = a/a ~ 7L 
is of the size of the discreteness scale, certainly an extreme regime in cosmology. Yet 
another intuitive way of expressing this regime is via densities: holonomy corrections 
are large when the matter density is of the order of the quantum-gravity density. By 
the classical Friedmann equation, this happens when 

We introduce the parameter 5hoi := p/Pqg = SnGv'^^^p/S in order to quantify holon- 
omy corrections. These are small when (5hoi *^ 1- 

The discreteness of loop quantum gravity manifests itself in different ways, some 
of which require more details to be derived. In addition to holonomy corrections, 
the most important one arises when one considers the inverse of the elementary 
lattice areas. Classically, the areas correspond to the densitized triad E^, which 
determines the spatial metric qab via -Ef-Ef = g"* det q and is canonically conjugate to 
the connection A^. The inverse of E'^ or its determinant appears in the Hamiltonian 
constraint of gravity as well as in all the usual matter Hamiltonians, especially in 
kinetic terms, and is thus crucial for the dynamics. 

Upon quantization, however, the densitized triad is represented in terms of the 
spin labels that also determine the lattice areas, and those labels can take the value 
zero. No densely defined inverses of the area operators exist, and therefore there 
is no direct way to quantize inverse triads or inverse volumes as they appear in 
Hamiltonians. However, as with holonomies replacing connection components, there 
is an indirect way of constructing well-defined inverse- volume operators, which imply 
further quantum corrections. 

2.2 Derivation of inverse- volume corrections 

The quantization of different kinds of inverse volumes or the co-triad ej^, obtained 



from the inverse of E^, begins with Poisson identities such as |l21| , P2 

Al / <l^x^\detE\ = 2Ti^Ge'^ha,c^=^=^ sgn(det E) = Ar^^Ge^ , (2.3) 
V J ) ^^\aeiE\ 

stemming from the basic Poisson brackets {^^(a;), Ej{y)} = SnG-ySjSaS^x, y). On the 
right-hand side of eq. ( p.3|) , there is an inverse of the determinant of E^, but on the 
left-hand side no such inverse is required. Classically, the inverse arises from deriva- 
tives contained in the Poisson bracket, but after quantization the Poisson bracket 
is replaced by a commutator and no derivative or inverse appears. In this way, one 
obtains well-defined operators for the inverse volume, implementing an automatic 
ultraviolet cutoff at small length scales. 
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The volume J d^x^/ldetEl of some region, containing the point v where we 
want to evaluate the co-triad, is quantized by well-defined volume operators, and the 
connection can be represented in terms of holonomies. For holonomies with edges of 
comoving length i^, we can write 

tr(r^/i.,e[/^-L K]) ~ hhioi^.je'' • (2.4) 

Here, r-' = ia^ /2 are Pauli matrices, /iv,e is a holonomy starting at a lattice vertex 
V in some direction e, and V^ is the volume of some region around v, with V^ its 
quantization. As long as v is included in the region integrated over to obtain the 
volume, it does not matter how far the region extends beyond v. One could even use 
the volume of the whole space. 

To quantize, loop quantum gravity provides the holonomy-fiux representation of 
the basic operators /iv,e (holonomies along edges e) and Fs = Jg d^y E^Ua, fiuxes 
of the densitized triad through surfaces 5* with co- normal Ua- These variables 
are SU (2)- valued, but one can devise a regular lattice for a simple implementa- 
tion of inhomogeneity, setting edges with tangent vectors e" = 5", / = 1,2,3 
in Cartesian coordinates. Then, holonomies are given by /iv.e^ = exp(£o''7c) = 
cos(£oc/2) -|- 2r/ sin(£oc/2) G SU(2), where c is the connection evaluated somewhere 
on the edge. All connection-dependent matrix elements can thus be expressed in 
the complete set of functions t] := exp(i£oc/2) G U(l), and the fiux through an ele- 
mentary lattice site in a nearly isotropic geometry is simply F = i^p with \p\ ~ a^, 
and p carrying a sign amounting to the orientation of space. Isotropy thus allows a 
reduction from SU(2) to U(l), with certain technical simplifications. 

For a nearly isotropic configuration, we assign a copy of the isotropic quantum 
theory to every (oriented) link / of a regular graph, making the theory inhomoge- 
neous. By this step we certainly do not reach the full theory of loop quantum gravity, 
which is based on irregular graphs with SU(2)-theories on its links. But we will be 
able to capture the main effects which have appeared in approximate considerations 
of loop quantum gravity with simpler graphs and reduced gauge groups. The basic 
operators are then a copy of fj^j and F^j for each lattice link with 

[57v,/, -Fv',j] = -47r74i'7v,/^7-j^v,v' , (2.5) 

if the edge of the holonomy and the surface of the flux intersect. 

When we insert holonomies for nearly isotropic connections in eq. ( p.4| ) and 
evaluate the trace, inverse-volume operators resulting from commutators have the 
form 

B.,i = 1^^:^ {vIjV.V.,1 - V.Avl,i) ■ (2.6) 

The volume at vertex v is obtained from components of the densitized triad, quan- 
tized by a flux operator F^j, with v an endpoint of the link J. If (/, /', /") denotes the 
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triple of independent links emanating from a given vertex, we can write the volume 



as K = \ \F^.iF^j'F^.i"\- Thus, 



B.,i = ^^^ (tiIi\J\KiF^J'Ki''\v^,i - V.j\l\KiKi^^\vl,i] ■ (2-7) 

As in the general representation, the basic operators F^j and f}v,7 satisfy the commu- 
tator identity ( |2.5|) while fj^j commutes with F^j/ and F^j". Moreover, fj^j satisfies 
the reality condition fj^jfjl j = 1. It turns out that these identities are sufficient to 
derive the form of inverse-triad corrections in a semiclassical expansion, irrespective 
of what state is used beyond general requirements of semiclassicality. 

We consider the two operators 'fi^,i\F^j\^^'^fil j and % j\F^j\^^'^fi^j separately.^ 
They can be simplified using the basic commutators ( p.5|) and [fjl j, F^j] = Air^yipifjl j. 



We can thus reorder terms so as to bring f/^ / right next to f/J j, and then cancel them 
using the reality condition. Reordering according to fj^jF^j = {F^,i — 47r7£pj)i7v^/ 
and fjl jF^j = {F^j + 4:nj£p{)i]l j leads to 



Vv,I 



\F.,i\'^^vli = \F.,i-'^^l4i\'^^, filj\F,,i\'/%j=\F,j + A7T^il,\'/\ (2.8) 



In the classical limit ^ — )■ 0, these expressions in B^j result in a derivative by F^j, 
as required by the Poisson bracket relationship ( p.3| ). For inverse- volume effects we 
are interested in the leading quantum corrections with h ^ 0, which arise in different 
forms. First, because the operator is nonlinear in the basic ones F^j and fj^j, classical 
expressions will be corrected by terms involving the moments of a state: As always in 
quantum physics, the expectation value (-Bv,/) does not have the classically expected 
relationship with expectation values of the basic operators. We can compute these 
corrections by following the principles of canonical effective dynamics, substituting 
(Fvj) + {Fvj — (Fvj)) for Fv,/ and performing a formal expansion by F^j — {Fv,i)'- 



(2.9) 



(2.10) 
(This expansion can be made well-defined and analyzed in the context of Poisson 



fl t\F rP/2^"^ 


= F^j-4n^il,\'/' 






-t p 1/2- 


= F.,/ + 4vr74ir/' 




^(^w)+4.74.^/^f:T ^^--^^-^^^ 



geometry of algebraic state spaces [^ g4|, P5[|.) 



^The factors of \J\F^j,\ and y |^v,/"| quantize positive powers of the densitized triad and do 
not give rise to inverse- volume corrections. 



All terms in F^j — (F^j) will either vanish upon taking an expectation value 
(A; = 1 in the expansion) or give rise to moments of the quantum state used to 
compute the expectation value (/c > 2, with fluctuations and correlations arising for 
k = 2). The precise values of the moments and their dynamics depend on the state 
used and in fact encode the state dependence of the theory, but for a semiclassical 
state they satisfy, by definition, the hierarchy {{F^j — (F^j))'') ~ /i^/^. This notion 
of semiclassicality is a very general one; it does not require us to choose a particular 
shape of the state, such as a Gaussian. 

The moment terms imply an important form of quantum corrections in the con- 
text of quantum back-reaction. Such corrections arise from different sources in the 
Hamiltonians, which will all have to be combined and analyzed. We will not enter 
such an analysis here, but rather note that even if we disregard quantum back- 
reaction, quantum corrections do remain: we have 

I f^t \p 11/2-. -. IP |l/2'.t \ 

47r74i 

where the dots indicate moment terms dropped. This expression includes inverse- 
volume corrections, computed for general semiclassical states. It matches with ex- 
pressions derived directly from triad eigenstates p6|, p7|, which are not semiclassical 



but, as proven here, provide reliable information about inverse- volume corrections. 
More general semiclassical states do not introduce additional dependence of inverse- 
volume corrections on connection components or curvature, they just introduce mo- 
ment terms which contribute to quantum back-reaction. (Such an extra dependence 
may arise from non-Abelian properties of the theory p8[ , which are not strong for 
perturbative inhomogeneities.) 

2.3 Correction functions 

Corrections to classical Hamiltonians in which inverse triad components appear can 
be captured by introducing correction functions such as 

a{a) := ^—^^ L{a) , (2.12) 

obtained by identifying 

(F.,,) = L2(a) (2.13) 

with the discreteness scale (depending on the scale factor in the presence of lattice 
refinement ||^). The multiplication of inverse- volume corrections by Lia) ensures 
that a[a) ~ 1 in the classical limit, but strong corrections can arise for small L. Our 
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derivations apply to small deviations from the classical value, for which consistent 
implementations in the dynamics are available. We can thus expand 

a{a) = 1 + ao^pi + ■ ■ ■ , (2.14) 

with 6pi := {£pi/L)"^ for m = A in the above derivation, and the dots indicating 
powers higher than m. 

For (-Fv,/) ^ ^pi inverse-volume corrections become very small, but they are 
significant if (-FV,/) is about as large as a Planck area or smaller. Bringing in our 
discreteness scale, leading inverse-volume corrections can be expressed in terms of 
the quantity 

(using m = 4 from now on). If L or f is constant, 6pi is constant and inverse- volume 
corrections merely amount to rescaling some expressions in Hamiltonians. More 
generally, however, the dynamical nature of a discrete state suggests that L and v 
change in time or, in cosmology, with respect to the scale factor a. We parameterize 
this dependence as 

Spioca~^ (2.16) 

with 0" > 0; see ^ for a discussion of possible values of a and its relation to quanti- 
zation parameters. 

In order to compare inverse-volume with holonomy corrections, we write 

The second equality shows that inverse- volume corrections are considerable and of the 
order one when the quantum-gravity density is close to the Planck density. Inverse- 
volume corrections thus behave very differently from what is normally expected for 
quantum gravity, where the Planck density is often presupposed as the quantum- 
gravity scale. In loop quantum gravity, this scale must be sufficiently small compared 
to the Planck density in order to be consistent with inverse-volume corrections. 
The last expression in eq. ( |2.17| ) is useful in order to compare holonomy with 



inverse-volume corrections. Inverse-volume corrections are usually suppressed by 
a factor of p/ppi, as expected for quantum-gravity effects, but there is an extra 
factor of 5j^qi. For small densities, holonomy corrections are small, but inverse- volume 
corrections may still be large because they are magnified by the inverse of (5hoi- As 
the energy density decreases in an expanding universe, holonomy corrections fall 
to small values, and in this way begin to magnify inverse-volume corrections. For 
instance, in an inflationary regime with a typical energy scale of p ~ 10~^'^ppi, we 
can use ( p.l7| ) to write (5hoi ~ 10^^/-\/5pi. Having small holonomy corrections of size 



5hoi < 10 then requires inverse-volume correction larger than 6pi > 10 . This 
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interplay of holonomy and inverse-volume corrections makes loop quantum gravity 
testable because it leaves only a finite window for consistent parameter values, rather 
than just providing Planckian upper bounds. It also shows that inverse-volume 
corrections become dominant for sufficiently small densities, as they are realized 
even in high-energy scenarios of inflation. 

In this context, it is worthwhile to comment on a comparison of the correc- 
tions derived here, assuming a nearly isotropic but explicitly inhomogeneous discrete 
state, with their form in pure minisuperspace quantizations. In inverse-volume as 
well as holonomy corrections, we referred to elementary building blocks of a discrete 
state, the plaquette areas in inverse-volume corrections and edge lengths in holon- 
omy corrections. A pure minisuperspace quantization would primarily make use of 
macroscopic parameters such as the volume of some region (or the scale factor) . The 
number of discrete blocks, such as M introduced above, is not available, and thus it 
is more difficult to refer to local microscopic quantities such as F^j. 

For curvature or the Hubble parameter, local quantities are easier to introduce 
and to use in holonomy corrections, but inverse-volume expressions must refer to 
quantities of size, which cannot be expressed microscopically in a pure minisuper- 
space context. As a consequence, inverse- volume corrections have often been misrep- 
resented in loop quantum cosmology. Without referring to A/", as it is introduced in 
the lattice-refinement formulation of loop quantum cosmology, one can only use the 
macroscopic volume of some region instead of the microscopic F^j.^ Inverse- volume 
corrections become smaller for larger -F^ /, and thus substituting this quantity by a 
macroscopic size suppresses the corrections. Any such suppression is merely an arti- 
fact of using the wrong expressions for the corrections based solely on minisuperspace 
considerations. Using a macroscopic volume also makes the corrections dependent on 
the size of the chosen region, which is another artificial dependence on extra param- 
eters; because of this, LQC inverse-volume corrections have been often interpreted 
as problematic or even unphysical. The derivation shown here solves these problems; 
see also the following subsection. 

As already seen, inverse-volume corrections show unexpected properties in terms 
of their dependence on the density, and regimes in which they are strong. Another 
unexpected property is seen in their infiuence on spacetime structure, with important 
consequences for cosmological perturbation theory. Inverse-volume corrections are 
not just of higher-curvature type in an effective action, but they deform the usual 
gauge algebra of generally covariant systems, generating spacetime diffeomorphisms. 
This deformation, as discussed in more detail in the following calculations, leads to 



^As mentioned earlier, in the inhomogeneous theory we can use the full volume or the size 
of any region in inverse-volume corrections because most plaquette contributions, which do not 
intersect the edge of the holonomies used, drop out. In homogeneous models, on the other hand, 
all plaquettes are equivalent and correspond to the same degree of freedom. The choices must thus 
be specified carefully in order to avoid minisuperspace artifacts. 
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characteristic cosmological effects. In a conceptual context, moreover, it allows us 
to distinguish inverse-volume corrections from the other types encountered in loop 
quantum gravity: holonomy corrections and quantum back-reaction. 

A closer look at the algebra of constraints generating the gauge transformations 
reveals that deformations of the algebra introduced by inverse-volume corrections 
cannot be undone by including holonomy corrections or quantum back-reaction [^. 
Holonomy corrections imply higher-order terms in the constraints depending on the 
connection nonpolynomially, or at least on the background connection if an expansion 
by inhomogeneities is done. No such terms arise for inverse-volume corrections, and 
no cancellation is possible. Quantum back-reaction, on the other hand, comes from 
terms including moments of a state, as alluded to in our derivation of inverse- volume 
corrections. The dependence on the moments remains if one computes the constraint 
algebra, in such a way that corrections from quantum back-reaction cannot cancel 
deformations implied by inverse-volume corrections, either. Since the characteristic 
effects analyzed here are a consequence of nontrivial deformations of the algebra, we 
can safely conclude that including only inverse- volume corrections does give a reliable 
picture, because they cannot be cancelled by the other, more complicated corrections. 
Of course, it remains of interest to study the inclusion of other effects such as the 
curvature of the universe, and the simultaneous competition between inverse- volume 
and other quantum corrections in a more complete dynamical analysis. 

2.4 Consistency 

Most of the properties and consequences of inverse- volume corrections are unexpected 
and unfamiliar. It is then perhaps not surprising that there are at least four main 
objections to the physical significance of effective LQC dynamics with inverse- volume 
corrections, which are popularly encountered in the literature and in scientific de- 
bates. It is claimed that (i) these corrections are ill-defined in a pure minisuperspace 
context and a flat universe, (ii) no rigorous derivation in the more involved inhomo- 
geneous context (taking into account lattice refinement) has been provided so far, 
(iii) even if a derivation were possible, the inflationary energy scale would be too low 
for volume/curvature corrections to be sizable, and (iv) even setting aside the issue 
of their size, the analysis would remain incomplete because we do not know how 
these corrections compete with holonomy modifications of the dynamics. As an ex- 
ample for the claimed incompleteness of correction functions used, the independence 
of inverse-volume corrections of the connection or curvature has been criticized as 
physically unjustified. 

We had already partially answered some of these objections elsewhere p. First, 
let us summarize the main arguments advanced there: 

(i) In a realistic cosmological scenario, there is no conformal invariance of the scale 
factor and the correct way to implement the quantum dynamics is to consider 
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the natural cell subdivision of space and how these cells evolve in time: this is 
the lattice refinement picture. In this perspective, interpretational difficulties 
regarding quantum corrections appear to be just an artifact of the idealized 
homogeneous and isotropic setting of pure minisuperspace models. 

(ii) Although a rigorous derivation is desirable, the motivations of lattice refine- 
ment are natural in the perspective of the full quantum theory and there is no 
conceptual obstacle in relaxing the parametrization obtained in a pure minisu- 
perspace.^ Moreover, one cannot simply suppress inverse-volume corrections 
by a regularization procedure, as occasionally suggested by taking the limit 
of Vo — ;■ C)0 in cases where these corrections are Vo-dependent. Inverse- volume 
corrections do appear in the full quantum theory and play an important role for 
well-defined Hamiltonians. If they disappeared by a regularization procedure 
in minisuperspace models, one should explain why they are absent in a cosmo- 
logical setting but not otherwise. Furthermore, there is tension between the 
requirement of closure of the inhomogeneous constraint algebra and the min- 
isuperspace parametrization [^, which demands clarifications; although the 
lattice parametrization is so far implemented semi-heuristically in calculations 
of effective constraint algebras, it does accommodate anomaly cancellation. 

(iii) Since the gauge symmetry of the model is deformed by quantum corrections, 
the very structure of spacetime is modiffed locally but everywhere; thus, one 
expects effects larger than in traditional scenarios of standard general relativ- 
ity with higher-order curvature terms. In [^ we found qualitative theoretical 
estimates of these effects which are several orders of magnitude larger than min- 
isuperspace estimates (and, interestingly, rather close to experimental bounds 
1^ in terms of orders of magnitude). However, the lack of control over the puta- 
tive quantum gravity characteristic scale (hidden in the quantum corrections) 
makes it difficult to assess its importance within inffation. 

(iv) We argued that other quantum corrections would not cancel inverse-volume 
effect because of the radically different way in which they affect the dynamics. 
Of course, the issue of comparing inverse- volume and holonomy corrections 
remains of interest for the community, but one does not expect that miraculous 
cancellations happen between the two. 



** Sometimes, the argum.ent is advanced that the minisuperspace parametrization (in particular, 
the so-called improved dynamics) is the only one producing a constant critical density and a robust 
bounce picture. This argument is invalid for two reasons. On one hand, even the improved dynamics 
parametrization does not give a constant critical density unless quantization ambiguities are tuned 
to certain specific values [[7| ; the time-dependent modification comes from inverse- volume corrections 
in the gravitational sector, which are nonzero in general. On the other hand, within the lattice 
parametrization a constant critical density, if desired, can be obtained, indeed. 



-13 



The results of the present section serve to further address the above objections 
and provide final clarifications for several of them. For the first time, we have embed- 
ded inverse- volume corrections in inhomogeneous models, using the lattice refinement 
picture and working at the kinematical level, thus giving fresh insight to these issues. 
In particular: 

i)-(ii) When the phase space volume is associated with an individual homogeneous 
cell rather than a fiducial volume (as done in pure minisuperspace) , the lat- 
tice parametrization emerges naturally and a quantum-gravity scale replaces 
unphysical quantities in inverse-volume corrections. Correction functions are 
completely independent of comoving volumes such as Vo and there is no regu- 
larization needed to make them disappear. 

(iii) Surprisingly, the magnitude of these corrections can be argued to be large at 
mesoscopic scales, even when densities are far away from Planckian values such 
as during inflation. In cosmological models, quantum corrections are relevant 
not just near a bounce at Planckian density. 

(iv) The basic noncancellation between inverse-volume effects and other, presently 
uncontrolled quantum corrections is reiterated with novel arguments. Modifi- 
cations of the classical constraint algebra by inverse- volume corrections cannot 
cancel with those from holonomy corrections, nor with terms from quantum 
back-reaction. Holonomy corrections provide an additional connection depen- 
dence of almost-periodic type in the constraints, while inverse-volume correc- 
tions as shown here have only weak connection dependence. Inverse-triad cor- 
rections are also independent of moments of a state, as they would determine 
quantum back-reaction. The structure of the Poisson algebra on the quantum 
phase space, including expectation values and moments, shows that neither the 
connection-dependent terms of the form of holonomy corrections nor moment 
terms describing quantum back-reaction can cancel the terms of inverse- volume 
corrections. If the constraint algebra is modified by inverse- volume corrections, 
it must remain modified when all corrections are included. Thus, also the pres- 
ence of effects larger than usually expected in quantum gravity is general. 

To summarize, loop quantum cosmology implies the presence of inverse-volume cor- 
rections in its cosmological perturbation equations. In their general parametrization, 
the corrections depend only on triad variables simply because they depend on a quan- 
tum scale whose dynamical nature is encoded by the background scale factor. This 
conclusion is a result of the derivations presented here, not an assumption. Also, 
their power-law form as a function of the scale factor is suggested by very general 
semiclassical considerations which do not further restrict the class of states. 
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3. Inflationary observables 

With a consistent implementation of inverse-volume corrections at hand, a com- 
plete set of cosmological perturbation equations follows. These equations have been 
derived elsewhere |]T1-|07 starting with a constraint analysis. Here we continue to 
prepare these equations for a convenient cosmological investigation, which we then 
exploit to find observational bounds on some parameters. 

The slow roll approximation is assumed precisely for the same reasons as in 
standard inflation: it is an Ansatz, it is the definition itself of inflation. In LQC 
there exists also a super-inflationary regime where the universe super-accelerates 
because of purely geometric effects, and the background attractor is not de Sitter. 
However, in that regime the constraint algebra has not been shown to close, and we 
have no rigorous control over the ensuing physics. Here, we simply assume that (i) 
inflation takes place thanks to a scalar field slowly rolling down its potential, and (ii) 
that this happens completely in the large-volume regime, where quantum corrections 
are small (and the algebra closes). 

3.1 Hubble slow-roll tower 

The slow-roll parameters as functions of the Hubble rate are defined starting from the 
background equations of motion, which also determine the coefficients of the linear 
perturbation equations. In the presence of inverse-volume corrections, the effective 
Friedmann and Klein-Gordon equations read 



«^4" 



^+PV(^) 



(3.1) 



and 

f" + 2-H n - ^^j if' + upV^ = , (3.2) 

respectively, where primes denote derivatives with respect to conformal time r : = 
J dt/a, 7i := a' /a = aH, k"^ = SnG, p = a"^, G is Newton's constant, and ^9 is a real 
scalar field with potential V{ip). Following section |^, the LQC correction functions 
are of the form 

a = 1 + ao^pi, (3-3) 

z/ = 1 + i/o<^Pi, (3.4) 

where ao and uq are constants and 

Spi oc a-" (3.5) 



is a quantum correction (2.16) whose time dependence is modelled as a power of the 



scale factor (here cr > is another constant). The proportionality factor will never 
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enter the analysis explicitly but, in the derivation of the perturbation equations, it is 
assumed that Spi < 1. Consistently, throughout the paper we use the equality symbol 
= for expressions valid up to 0{6pi) terms, while we employ ~ for relations where the 
slow-roll approximation has been used. The latter holds when the following slow-roll 
parameters are small: 



e := 1 
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1 + 



"0 + t'o ( g - 1 



a 
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The conformal-time derivatives of e and t] are 

e' = 2He{e - r]) - aUiSpi , 



where 



e:=ao[- + 2e-r] 
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1 e 



(3.6) 
(3.7) 

(3.8) 
(3.9) 

(3.10) 



The inflationary spectra were computed in Pj. The scalar power spectrum is 

(3.11) 
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Equation ( p.ll ) is evaluated at the time k = H when the perturbation with comoving 
wavenumber k crosses the Hubble horizon. Using the fact that 

6'p^ = ~aH6pr (3.13) 

and ' ~ Tid/dlnk, the scalar spectral index Ug — 1 := dlnPs/dlnfc reads 

n, - I = 2r] - Ae + a^njpu (3.14) 

where 

7n, := - - ao ( 1 ) - 7s = "0 - 2z/o H — - , (3.15) 

eve/ a+1 

while the running Og := dng/dlnA; is 

as = 2{5er] - Ae^ - e) + ^(4e - a^nJSpi . (3.16) 
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This shows that, for a = 0(1), the running can be as large as 6p\. In this case, the 
terms higher than the running can give rise to the contribution of the order of Sp\. 



In section ^?1| we shall address this issue properly. 
The tensor power spectrum is 

Vt = —{l+-ft6pi), 



7t 



a- 1 



«0' 



(3.17) 



while its index rit := dlnPt/dlnfc and running at := drit/dlnk are, respectively, 

fit = -2e - cr'jtSpi , (3.18) 

and 

at = -4e(e - r]) + a{2i + a^t)Spi . (3. 19) 

The tensor-to-scalar ratio r := Vt/Vg combines with the tensor index into a consis- 
tency relation: 

r = 16e [1 + (7t - 7.)5pi] 
= -8{n, + [nt(7t - 7s) + a7t]5pi} • (3.20) 

When 6pi = 0, all the above formulas agree with the standard classical scenario [^ . 

3.2 Potential slow-roll tower 

To constrain the inflationary potential against observations, it is convenient to recast 
the cosmological observables in terms of the tower of slow-roll parameters written as 



functions of V and its derivatives. From eqs. ( p.lj) and ( p.6[) we have 
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The first three elements of the tower are 
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On using eqs. ( 3.23|) - (|3.25|) , we have the following technical expressions 
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which give the inversion formulas 
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We can now rewrite the cosmological observables. The scalar index ( 3.14| ) and 
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its running ( p.l6 ) become 

Hs-l 
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The tensor index ( p.l8|) and its running ( |3.19|) are 



Ut = -2ey -c„,5pi, 

at = -4ev (2ev. - riy) + c^t^pi 
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Finally, the tensor-to-scalar ratio ( p.20 ) is 

r = 16ev + CrSpi 



where 



8[3ao(3 + 5a + Ga^) - z/o(t(6 + 11a)] IGaao 



9(a+l) 



(3.37) 



(3.38) 



4. Power spectra and cosmic variance 

In this section, we cast the power spectra as nonperturbative functions of the wavenum- 
ber k and a pivot scale ko (section ^3]) . The parameter space of the numerical analysis 
is introduced in section [4.2| , while a theoretical prior on the size of the quantum cor- 



rection is discussed in section [4.3| . An important question to address is whether a 
possible LQC signal at large scales would be stronger than cosmic variance, which is 
the dominant effect at low multipoles. This issue is considered in section [4.4| , where 
a positive answer is given for a certain range in the parameter space. 

4.1 Power spectra and pivot scales 



Because of eq. (|3.13| ), terms higher than the runnings Og and at can give rise to a 
nonnegligible contribution to the power spectra Vs{k) and Vt{k). Let us expand the 
scalar power spectrum to all orders in the perturbation wavenumber about a pivot 
scale ko'- 



\nV,{k)=\nV,{ko) + [n,{kc 



Ti , "s(A;o) 2 

l\x -\ X 






m=3 
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(4.1) 



where x := ln(A;/fco), and as := d'"~^as/(dln /c)™"^. When 0{ev) and 0{r]v) terms 
are ignored, c„^ ~ /s in eq. ( p.29|) , while the dominant contribution to the scalar 
running can be estimated as 
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k=ko 



Similarly, we can derive the m-th order terms as as 

a^)(A;o)^(-l)'"a"-Vs5pi(A;o)- (4-3) 

In this case, the last term in eq. ( W^ converges to the exponential series, 
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Thus, the scalar power spectrum (4J.) can be written in the form 



V,{k) = Vsiko) exp <; [n^iko) - l]x + ^^x' 



+fsSpi{ko 



x\l- lax] + -(e-"^' - 1) 
2 / cr 



(4.5) 



We stress that the power spectrum ( [4.5|) is nonperturbative in the wavenumber k, 
while ordinary inflationary analyses take a truncation of ( [4.11) . The expression (|4.5| ) 
is valid for any value of a and of the pivot wavenumber, provided the latter lies within 
the observational range of the experiment. Note that ko is not fixed observationally, 
except from the fact that we can choose any value on the scales relevant to CMB 
(with the multipoles i ranging in the region 2 < i < 1000). The CMB multipoles 
are related to the wavenumber k by the approximate relation 

k^ IQ-^hi Mpc'S (4.6) 

where we take the value /i = 0.7 for the reduced Hubble constant. The default pivot 
value of CMBFAST [|0l and camb |^ codes is fco = 0.05 Mpc^^ (4 ~ 730). For 



the WMAP pivot scale ko = 0.002 Mpc"^ {£ ~ 29) [|3], H], the maximum value 
of X relevant to the CMB anisotropics is Xmax ~ 3.6. Intermediate values of ko are 



also possible, for instance ko = 0.01 Mpc~ |3^. In general, the constraints on the 
parameter space, and in particular the likelihood contours, depend (even strongly) 
on the choice of the pivot scale ||36[, and it is interesting to compare results with 
different ko also in LQC. 

The fact that we can resum the whole series is of utmost importance for the 
consistency of the numerical analysis. In standard inflation, higher-order terms do 
not contribute to the power spectrum because they are higher-order in the slow-roll 
parameters. Then, one can truncate eq. ( [4.1|) to the first three terms and ignore 
the others. Here, on the other hand, all the terms ( [4.3|) are linear in 6pi and they 
contribute equally if the parameter a is large enough, a > 1. This fact might naively 
suggest that small values of a are preferred for a consistent analysis of a quasi-scale- 
invariant spectrum P). In that case, one would have to impose conditions such as 
|[ns(A;o) — l]a;| ^ |[tts(fco)/2]a;^|, which depend on the pivot scale ko- 

For cr > 1, however, different choices of ko would result in different convergence 
properties of the Taylor expansion of Vs- The point is that Sp\{k) changes fast for 
cr > 1 and the running of the spectral index can be sizable; dropping higher-order 
terms would eventually lead to inconsistent results. On the other hand, eq. (|4.5| ) 
does not suffer from any of the above limitations and problems, and it will be the 



basis of our analysis, where ns(fco) and as{ko) are given by eqs. ( |3.29|) and (|3.30|) 



The last term in eq. ( [4.5| ), usually negative, tends to compensate the large positive 
running, thus providing a natural scale-invariance mechanism without putting any 
numerical priors. 
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Assuming that c„j ~ ft, same considerations hold for the tensor spectrum, which 
can be written as 

«t(fco) 2 



Vt{k) = Vt{kQ) exp <^ nt{kQ) x H — x 



+ftSpi{ko) 



x{l- -ax] + -(e-"^^ - 1) 



(4.7) 



where nt(/i;o) and at(^o) are given by eqs. ( |3.32| ) and ( p.33| ), respectively. Finally, the 
tensor-to-scalar ratio is given by eqs. ( |3.37] ) and ( p. 381 ), with the slow-roll parameters 
evaluated at the pivot scale k = k^. 

4.2 Parameter space 

The CMB likelihood analysis can be carried out by using eqs. ( [4.5|) , ([4.7|), and ( |3.37| ). 
Let us take the power-law potential [p7| 

V{ip) = Voip''. (4.8) 

In this case, it follows that (/cq dependence implicit) 

- - "' „.^H^^i^.., ,;^ 4(»-i)("-^) 4. (4.9) 



This allows us to reduce the slow-roll parameters to one (i.e., e^). 
For the exponential potential [^ 

\/(<^) = K)e-''^^ (4.10) 

the relation between the slow-roll parameters is given by 

ev = y, W = 2e^, e' = 4e^ , (4.11) 

which are again written in terms of the single parameter e^. 

Between the model parameters vq and ao we can also impose the following rela- 
tion p, valid for cr 7^ 3: 

[a + <6)[a - 3) 
Introducing the variable 

<5(A;o) := ao5pi(M , (4-13) 

we can write /s(5pi(A;o) and /t(5pi(A;o) in the form 

,,,,. a(8a3-8a2-93(T+18)^,, , ^ ^ ,, ^ 2^2 

/s5pi(fco) = ^j7 ^w ^,w ^./ '^(fco) , /t5pi(A;o) = ^— 5(A;o) . (4.14) 

2((T — 3j(cr -|- lj(cr -|- oj a + 1 

For cr = 3 one has Oq = identically, in which case eq. ( |4.13|) is replaced by diko) := 
i^oSpiiko). 

To summarize, using the relation ( [1.12| ), all the other observables can be written 
in terms of (5(fco) and ey(A;o)- Hence, for given a and ko, one can perform the CMB 
likelihood analysis by varying the two parameters 6{ko) and ev{ko). 
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4.3 Theoretical upper bound on the quantum correction 

For the validity of the linear expansion of the correction functions ( p.3| ) and ( |3.4| )^ and 
all the perturbation formulas where the 0{6pi) truncation has been systematically 
implemented, we require that 6{k) = ao^Pi(^) < 1 for all wavenumbers relevant to 
the CMB anisotropics. Since 5pi oc a~'^, the quantity 6{k) appearing in inflationary 
observables is approximately given by 

S{k) = 6{ko)(^^y , (4.15) 

where we have used k = % aX Hubble exit with T-L/a ^ const. As k (x £, the 
same expression can be written in terms of the multipoles i. Since a > 0, one has 
6{k) > 6{kQ) for k < k^ and 5{k) < 6{kQ) for k > k^. This means that the larger the 
pivot scale /cq, the smaller the upper bound on 6{ko). 

Let us consider two pivot scales: (i) ko = 0.002 Mpc~^ (multipole io ~ 29) 
and (ii) ko = 0.05 Mpc~^ (multipole io ~ 730). Since the largest scale in CMB 
corresponds to the quadrupole £ = 2, the condition 6{k) < 1 at £ = 2 gives the 
following bounds (5max on the values of S{ko) with two pivot scales: 



i) 5^ax = 14.5-'^ 


(for ko = 0.002 Mpc~^) , 


(4.16) 


ii) <5max = 365-'^ 


(for ko = 0.05 Mpc"^) . 


(4.17) 



Values of (5max for some choices of a are reported in table |I|. The suppression of 
(^max for larger ko and a can be also seen in the power spectra ( [4. 5] ) and ([4.71). The 



term e""""^ = [ko/k)" can be very large for large ko'. for instance, if a = 6 and 
ko = 0.05 Mpc^^, one has e~°"^ ~ 10^^ at £ = 2. Then we require that S{ko) is 
suppressed as S{ko) < 10~^^. 

4.4 Cosmic variance 

At large scales, the failure of the ergodic theorem for the CMB multipole spectrum 
manifests itself in the phenomenon of cosmic variance, an intrinsic uncertainty on 
observations due to the small samples at low multipoles. For a power spectrum V{i), 
cosmic variance is given by 

Vs^Mi) = ^^V\i) . (4.18) 

A natural question, which is often overlooked in the literature of exotic cosmologies, is 
how effects coming from new physics compete with cosmic variance. In our particular 
case, we would like to find which values of a give rise to a theoretical upper bound 



^Using the relation (4.12), one sees that vq is of the same order as Uq, so a bound on 6 is 
sufficient. 
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a 


0.5 


1 


1.5 


2 


3 


6 


ko = 0.002 Mpc"^ 


'-'max 

6 


0.26 
0.27 


6.9 X 10-2 
3.5 X 10-2 


1.8 X 10-2 
1.7 X 10-3 


4.7 X 10-3 

6.8 X 10-5 


3.2 X 10-"^ 

4.3 X 10-^ 


1.0 X 10-^ 


ko = 0.05 Mpc"^ 


"max 

s 


5.2 X 10-2 
6.7 X 10-2 


2.7 X 10-3 
9.0 X 10-^ 


1.4 X 10-^ 
1.3 X 10-5 


7.5 X 10-6 
1.2 X 10-^ 


2.1 X 10-8 
2.7 X 10-" 


4.3 X 10-16 



Table 1: Theoretical priors on the upper bound of S (= (^max) and 95% CL upper limits 
of 5 constrained by observations for the potential V{(p) = V^^p^ with different values of a 
and for two pivot scales. The likelihood analysis has not been performed for a = 6 since 
the signal is below the cosmic variance threshold already when a = 2. For o" = 3, the 
parameter 5 = I'o'^Pi has been used. 



(^max of inverse-volume LQC quantum corrections larger than the error bars due to 
cosmic variance with respect to the classical spectrum. 

Consider the scalar spectrum Vs{i), eq. ( [4.5|) with k/k^ replaced by i/io- It is 
determined up to the normalization Psl^o), so that the region in the (£, Vs{i)/Vs{io)) 
plane affected by cosmic variance is roughly delimited by the two curves 



PsW ± VVarp,(£) 



n 



1± 



nw 



2£ + l / Vsiio) 



(4.19) 



where we take the classical spectrum as reference. 

The power spectrum ([4.5|) , together with the cosmic variance effect ( [4.19|) , is 
shown in figure |I| for n = 2 and the pivot scale io = 29. Ignoring the solid lines 
for the moment, the dashed curves correspond to 6 = (5max- The exponential term 
e"°"^ = [ko/kY in eq. ( |4.5| ) gives rise to an enhancement of the power spectra on large 
scales, as we see in the figures. ^ For a > 3, the growth of this term is so significant 
that 5{€) must be very much smaller than 1 for most of the scales observed in the 
CMB, in order to satisfy the bound 5{l = 2) < 1. LQC inverse- volume corrections 
are well within the cosmic variance region for cr > 2. However, already at a = 2 
quantum corrections strongly affect multipoles £ < 6. For a ~ 1, the spectrum 
is appreciably modified also at multipoles (i > 500. Changing the pivot scale to 
C.Q = 730, one sees that the quantum effect is generally greater than cosmic variance 
at sufficiently low multipoles (figure |]). The plots for n = 4 are very similar and we 
shall omit them. 



5. Likelihood analysis 

We carry out the CMB likelihood analysis for the power-law potential 
^This feature is similar to results for tensor modes fpll. 



as well as 
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Figure 1: Primordial scalar power spectrum Vs{i) for the case n = 2, with ev(fco) = 0.009 
and the pivot wavenumber fco = 0.002 Mpc~^, corresponding to £o = 29. The values of a 
are a = 2 (top panel) a = 1.5 (center panel), and a = 1 (bottom panel), while we choose 
three different values of 6{£o), as given in table |l]: (classical case, dotted lines), the 
observational upper bound from the numerical analysis (solid lines), and Jmax (a-priori 
upper bound, dashed lines). Shaded regions are affected by cosmic variance. 



the exponential potential ( |4.10|) . We run the Cosmological Monte Carlo (CosmoMC) 
code [^ with the data of WMAP 7yr |^ combined with large-scale structure (LSS) 

and Big Bang 



including BAO), HST |ET|, Supernovae type la (SN la) 
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Figure 2: Primordial scalar power spectrum Vs{i) for the case n = 2, with ev(fco) = 0.009 
and the pivot wavenumber ko = 0.05 Mpc~^, corresponding to io = 730. The values of a 
are a = 2 (top panel), o" = 1.5 (center panel), and a = 1 (bottom panel), while we choose 
three different values of S{£q), as given in table |l]: (classical case, dotted lines), the 
observational upper bound from the numerical analysis (solid lines), and 5max (a-priori 
upper bound, dashed lines). Shaded regions are affected by cosmic variance. 



Nucleosynthesis (BBN) |Q, by assuming a ACDM model. 

In the Monte Carlo routine we vary two inflationary parameters S{kQ) and ev{ko) 
as well as other cosmological parameters. Note that S^ko) and eviko) are constrained 
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at the chosen pivot scale ko. While the bound on 6 depends on ko (and it tends to 
be smaller for larger ko), that on (koYS^ko) does not [see eq. (|4.15|) ]. 



Under the conditions e^ ^ 1 and 6 '^ 1, the slow-roll parameter e^ is approxi- 
mately given by e^ ~ (fi:^/2)(y9'/'H)^. Then the number of e-foldings during inflation 
can be estimated as 

N:= rdrU^K r d^—l== (5.1) 

Jt Jiff ^y2ev{!f) 

where (/?£ is the field value at the end of inflation determined by the condition e^ ~ 
0(1). For the power-law potential (|4.8| ) one has (ft ~ n/v2K? and A^ ~ 'n,/(4ev)— n/4, 
which gives 

ey^—- , 45<A^<65. (5.2) 

AN + n ^ ' 

The typical values of A^ for the perturbations relevant to the CMB anisotropies are 
actually around 50 < A^ < 60, but we have taken the wider range above. The 
comparison of this estimate with the experimental range of ey will determine the 
acceptance or exclusion of an inflationary model for a given n. 

For the exponential potential ( [4.10|) the slow-roll parameter e^ is constant, which 
means that inflation does not end unless the shape of the potential changes after some 
epoch. In this case, we do not have constraints on ty coming from the information 
of the number of e-foldings in the observational range. 

5.1 Quadratic potential 

Let us study observational constraints in the case of the quadratic potential V{ip) = 

5.1.1 ko = 0.002 Mpc-^ 

We first take the pivot wavenumber ko = 0.002 Mpc"^ {io ~ 29) used by the WMAP 
team [^. In figure |], the 2D posterior distributions of the parameters S{ko) and 
ey{ko) are plotted for n = 2 and cr = 2, 1.5, 0.5. We have also run the code for other 
values of a such as 1 and 3. The observational upper bounds on S are given in table 
|I] for several different values of a. 

For a < 1, the exponential factor e''^^ does not change rapidly with smaller val- 
ues of /s,t5 so that the LQC effect on the power spectra would not be very significant 
even if 6{ko) was as large as ey{ko). As we see in figure |l] (solid curve), if cr = 0.5 
the LQC correction is constrained to be S{ko) < 0.27 (95% CL), which exceeds the 
theoretical prior Sma.^ = 0.26. Since S{ko) is as large as 1 in such cases, the validity 
of the approximation 6{ko) < ey{ko) to derive the power spectra is no longer reliable 
for a < 0.5. 

Looking at table |1|, when a = 1, the observational upper bound on 6{ko) becomes 
of the same order as 5max- Foi' o" < 1.5 the effect of the LQC correction to the power 
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Figure 3: 2-dimensional marginalized distribution for the quantum-gravity parameter 
5{ko) and the sfow-roll parameter ev(feo) with the pivot /cq = 0.002 Mpc~^ for n = 2, 
constrained by the joint data analysis of WMAP 7yr, LSS (including BAO), HST, SN la, 
and BBN. The values of a are a = 2 (top left panel), a = 1.5 (top right panel), and a = 0.5 
(bottom panel). The internal and external lines correspond to the 68% and 95% confidence 
levels, respectively. 



spectrum becomes important on large scales relative to cosmic variance. For smaller 
a the observational upper bound on 6{ko) = ao^Pi(^o) tends to be larger. When 
a = 1.5 the LQC correction is constrained to be S^ko) < 1.7 x 10~^ (95% CL), see 
figure ^ This is smaller than the theoretical prior 5max = 1.8 x 10~^ by one order of 
magnitude. 

The effect of cosmic variance is significant for a > 1.5. When o" = 3, the LQC 
correction is constrained to be S{ko) = t'o^Pi(^o) < 4.3 x 10^^ (95% CL). With respect 



to the prior 6^ 



3.3 X 10 '^, the observational bound is smaller by three orders of 



-28- 



magnitude. For cr > 3 the power spectra grow very sharply for low multipoles, so 
that the upper bounds on S{ko) become smaller. Numerically it is difficult to deal 
with such rapidly changing power spectra. 

For the sake of completeness, we should notice that the bounds plotted in Figs. |l| 
and 1^ include input from several datasets, but the cosmic variance belt comes only 
from the CMB. Therefore, the medium-scale part of the cosmic- variance plots might 
not give the full picture of the statistical limitations in this range, but we do not 
expect appreciable modifications from large-scale structure observations. 

For n = 2, the theoretically constrained region ( |5.2| ) corresponds to 0.008 < 
ey < 0.011. As we see in figure ^, for a > 0.5, the probability distributions of e^ are 
consistent with this range even in the presence of the LQC corrections. Hence, for 
the pivot wavenumber ko = 0.002 Mpc~ , the quadratic potential is compatible with 
observations as in standard cosmology. 



5.1.2 ko = 0.05 Mpc-i 



,x10 





Figure 4: 2-dimensional marginalized distribution for n = 2 with the pivot ko = 0.05 
Mpc^^. The values of a are a = 2 (left panel) and a = 1 (right panel). 

We proceed to the case of the pivot wavenumber ko = 0.05 Mpc~^ (io ~ 730). 
From eq. ( 4.17 ), the theoretical priors on (5max for given a are smaller than those 
corresponding to ko = 0.002 Mpc~^. In figure ^ we plot the 2-dimensional poste- 
rior distribution of S{ko) and ev{ko) for n = 2 and a = 2, 1. When a = 2, the 
observational upper limit is found to be S{ko) < 1.2 x 10~^ (95% CL), which is two 
orders of magnitude smaller than the bound 6{ko) < 6.8 x 10~^ obtained for the pivot 
ko = 0.002 Mpc~^. This comes from the fact that the choice of larger ko leads to 
more enhancement of power on large scales. 
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When cr = 1, we find the constraint 6{ko) < 9.0 x 10^"^ (95% CL), which is 
about 1/3 of the theoretical prior 5max = 2.7 x 10~^. From table |I|, we see that the 
observational limit of 6 ior a = 0.5 exceeds (5max- Hence, our combined slow-roll/5pi 
truncation is no longer trustable for a < 0.5, as it happens for kg = 0.002 Mpc~^. 

From figure ^ we find that the theoretically allowed range of ey (0.008 < ey < 
0.011) is consistent with its observational constraints. The different choice of k^ 
affects the upper bounds on ^(fco), but the basic property of the LQC effect on the 
power spectra is similar. The quadratic inflaton potential can be consistent with 
the combined observational constraints even in the presence of the LQC corrections, 
independent of the values of fco relevant to the CMB anisotropies. 

5.2 Quartic potential 
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Figure 5: 2-diinensional marginalized distribution for n = 4 in tlie two cases: (i) cr = 1 
and ko = 0.002 Mpc"^ (left panel) and (ii) a = 2 and ko = 0.05 Mpc"! (riglit panel). 



Let us proceed to the case of the quartic potential V{(p) = Voip'^. Numerically, 
we find that the observational upper bounds on 6{ko) for given a and ko are similar 
to those for the quadratic potential. In the top panel of figure |^ we show the 2- 
dimensional distribution for cr = 1 with the pivot wavenumber fco = 0.002 Mpc~^. 
The LQC correction is constrained to be S{ko) < 3.4 x 10~^ (95 % CL), which is 
similar to the bound S^ko) < 3.5 x 10~^ for n = 2 (see table jT]). The bottom panel 
of figure 1^ corresponds to the posterior distribution for a = 2 with ko = 0.05 Mpc~^, 
in which case 6{ko) < 1.1 x 10~^ (95 % CL). Since the LQC correction given in 
eq. ( |4.15|) does not depend on the values of n, the above property of ra-independence 
can be expected. For larger a and ko the upper bounds on 6{ko) tend to be smaller. 

From eq. ( ^.2|) , the values of e^ related with the CMB anisotropies fall in the 
range 0.015 < e^ < 0.022. For cr = 1 and ko = 0.002 Mpc~^, this range is outside 
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the 1(7 likelihood contour. In particular, for A^ < 58, this model is excluded at the 
95% confidence level. For o" = 2 the observationally allowed region of ev is slightly 
wider than that for a = 1. However, as we see in the lower panel of figure ^, this 
model is still under an observational pressure. Numerically we have confirmed that 
the bounds on ev{ko) are insensitive to the choice of k^. Hence the quartic potential 
is in tension with observations even in the presence of the LQC corrections. 

5.3 Exponential potentials 
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Figure 6: 2-dimensional marginalized distribution for the exponential potential V{}p) = 
Voe^'^-^'P in the two cases: (i) cr = 2 and ko = 0.002 Mpc"^ (left panel) and (ii) a = 1.5 and 
ko = 0.05 Mpc~^ (right panel). 

Finally, we study the case of exponential potentials. In figure ^ we show the 2- 
dimensional posterior distribution for (i) a = 2 and ko = 0.002 Mpc~^ and (ii) a = 1.5 
and ko = 0.05 Mpc~^. The observational upper limits on the LQC corrections for the 
cases (i) and (ii) are ^(fco) < 6.8 x 10~^ and (5(A;o) < 1-3 x 10~^, respectively, which 
are similar to those for n = 2 with same values of a. Hence, for given values of a and 
ko, the effect of the LQC corrections to the power spectra is practically independent 
of the choice of the inflaton potentials. 

On the other hand, the observational constraints on the slow-roll parameter 
depend on the potential. In figure ^ we find that the observationally allowed values 
of ev{ko) are in the range 0.005 < ev{ko) < 0.27 (95 % CL) for two different choices 
of kf). The maximum value of ev{ko) is larger than that for n = 2 and n = 4. Since 
inflation does not end for exponential potentials, one cannot estimate the range of 
the slow-roll parameter relevant to the CMB anisotropics. Hence one needs to find 
a mechanism of a graceful exit from infiation in order to address this issue properly. 
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6. Conclusions 

In the presence of the inverse- volume corrections in LQC, we have provided the ex- 
phcit forms of the scalar and tensor spectra convenient to confront inflationary models 
with observations. Even if the LQC corrections are small at the background level, 
they can significantly affect the runnings of spectral indices. We have consistently 
included the terms of order higher than the scalar /tensor runnings. Inverse-volume 
corrections generally lead to an enhancement of the power spectra at large scales. 

Using the recent observational data of WMAP 7yr combined with LSS, HST, 
SN la, and BBN, and analyzing them with techniques routinely used also in stan- 
dard inflation, we have placed constraints on the power-law potentials V{(p) = Vqv?" 
{n = 2,4) as well as the exponential potentials V{{p) = Voe"**^"^. The inflationary 
observables (the scalar and tensor power spectra Vs, Vt and the tensor-to-scalar ratio 
r) can be written in terms of the slow-roll parameter e^ = {V^/V)'^/{2k,'^) and the 
normalized LQC correction term 6. We have carried out a likelihood analysis by 
varying these two parameters as well as other cosmological parameters for two pivot 
wavenumbers fco (0.002 Mpc~^ and 0.05 Mpc~^). 

The observational upper bounds on 6{ko) tend to be smaller for larger values 
of kf). In table |I| we listed the observational upper limits on 6{ko) as well as the 
theoretical priors (Jmax for the quadratic potential V{(p) = Voip^ with a number of 
different values of the quantum gravity parameter a (which is related to (5 as 5 oc a""^). 
For larger a, we find that (5(A;o) needs to be suppressed more strongly to avoid the 
significant enhancement of the power spectra at large scales. When a < 0.5 the 
observational upper limits of (5(fco) exceed the theoretical prior (5max, which means 
that the expansion in terms of the inverse-volume corrections can be trustable for 
o^>0.5. 

As we see in Figs. ^^ and in table |I], the observational upper bounds on 5(/co) 
for given k^ and a are practically independent of the choice of the inflaton potentials. 
This property comes from the fact that the LQC correction for the wavenumber k is 
approximately given by 5{k) = S{ko){ko/ky, which only depends on ko and a. On 
the other hand the constraints on the slow-roll parameter e^ are different depending 
on the choice of the inflaton potentials. We have found that the quadratic potential 
is consistent with the current observational data even in the presence of the LQC 
corrections, but the quartic potential is under an observational pressure. For the 
exponential potentials the larger values of e^ are favored compared to the power- 
law potentials. However, the exponential potentials are not regarded as a realistic 
scenario unless there is a graceful exit from inflation. 

The exponential term e~'^^ = {ko/k)" in eq. ( [4.5|) is responsible for the enhance- 
ment of the power spectrum at large scales. This feature is characteristic of the model 
and is not reproduced by other sources. For instance, non-commutative geometry 
or string corrections P3 predict a suppression, rather than an enhancement, of the 
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spectra. Also some old papers on LQC advertized a suppression of power (e.g., the 
second reference of ||l3)! but the quantum corrections were not under full control at 



the level of perturbation theory; the present results supersede those early discussions. 
Moreover, the signatures of the LQC spectra cannot be mimicked by any standard 
scalar potential. The enhancement of power is due to a scale-dependent correction 
in the spectral amplitudes, while exotic potentials would not affect the perturbed 
dynamical equations (compare with the LQC Mukhanov equations |^). In particu- 
lar, the consistency relation (|3.2CI|) is notably different with respect to the standard 
relation r = — 8nt. Finally, also Wheeler-DeWitt quantum cosmology predicts an 



enhancement of power at large scales [45|, but the effect is qualitatively different 



from the structure of 6pi and its size is much smaller. The main reason is that it is 
governed by the energy scale of inflation, contrary to what happens in LQC. 

If we compare the observational upper bounds with the theoretical lower bounds 
discussed in section §, we can see that estimates of these parameters are separated by 
at most a few orders of magnitude, much less than is usually expected for quantum 
gravity. By accounting for fundamental spacetime effects that go beyond the usual 
higher-curvature corrections, quantum gravity thus comes much closer to falsifiability 
than often granted. It is of interest to see how the future high-precision observations 
such as PLANCK will constrain the LQC correction as well as the slow-roll param- 
eters. Even in the case where the quadratic potential were not favored in future 
observations, it would be possible that the small-field inflationary models be consis- 
tent with the data. For these general inflaton potentials, the effect of inverse- volume 
corrections on the CMB anisotropies should be similar to that studied in this paper. 
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